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ABSTRACT 


An aeroelastlc stability of three-degree 
flap-lag-torslon blade In forward flight Is 
examined. Quasisteady aerodynamics with a 
dynamic Inflow model Is used. The nonlinear 
time dependent periodic blade response Is 
calculated using an Iterative procedure based 
on Floquet theory. The periodic perturbation 
equations are solved for stability using 
Floquet transition matrix theory as well as 
constant coefficient approximation In the 
fixed reference frame. Results are presented 
for both stiff-1 nplane and soft-1 nplane blade 
configurations. The effects of several para- 
meters on blade stability are examined. 
Including structural coupling, pitch-flap and 
pitch-lag coupling, torsion stiffness, steady 
Inflow distribution, dynamic Inflow, blade 
response solution and constant coefficient 
approximation. 


NOTATIONS 


■ lift curve slope 

* matrix In first order 
equations In rotating system 

■ matrix In first order 
equations In fixed system 

• matrices In dynamic Inflow 
equations 

■ blade chord 
« blade section drag coef- 
ficient 

■ blade section lift coef- 
ficient 

■ bltde section moment coef- 
ficient 
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» differential aerodynamic 
coefficients wrt a 

• drag force of the helicopter 
» hinge offset divided by 

rotor radius 

3 equivalent drag area of 
helicopter 

3 nonlinear force vector in 
response equation 
3 moment Inertial of blade 
(flap) 

3 ratio of torsional Inertia 
to blade flap Inertia 
3 distance of hub from heli- 
copter c.g. 

• rotor drag force, positive 
rearward 

3 pitch-flap and pitch-lag 
couplings respectively 
3 stiffness matrices In response 
and stabllty equations 
3 coefficients In Drees model 
3 coefficient matrices In 
dynamic Inflow equation 
3 mass matrices In response 
and stability equations 
3 Aerodynamic rolling and 
pitching moments respectively 
3 aerodynamic flap, lag and 
pitch moments respectively 
3 number of blades 
3 matrix In fixed system 
defined by eqn. (29) 

3 transition matrix 

• rotor radius 

3 structural coupling para- 
meter 
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■ rotor thrust force 

» blade section normal 
and Inplane velocity 

* blade secti on resulta nt 
velocity, ^2 + ,j 2 

t p 

« forward speed 

* helicopter gross weight 

3 blade radial coordinate 

(nondlmenslonallzed wrt 
radius) 

« chordwise offset of blade 
aerodynamic center behind 
pitch axis 

= chordwise offset of blade 
center of gravity behind 
pitch axis 

■ vector consisting degrees of 
freedom In rotating system 

« state vector In rotating 
system 

■ rotor side force, positive 
towards advancing side 

« fuselage side force 

> blade section angle of attack 

» real part of kth charac- 
teristic exponent 

* angular deflections (flap, 
lag, torsion) 

= precone angle 

■ blade lock number, paclr/I b 

* collective, lateral cyclic, 
longitudinal cyclic pitch 
angles 

* rotor Inflow ratio 

» kth characteristic exponent 

* advance ratio, V cos a/ PR 

» rotating flap, lag and torsional 
frequencies 

* nonrotating flap, lag and tor- 
sional frequencies 

* air density 

* solidity ratio, N fa c/tR . 

» section Induced angle, tan * 1 u /U t 

* lateral tilt of shaft p 

* azimuth angle of the blade 

* Imaginary part of the kth exponent 

■ speed of rotation 

■ perturbation quantity 


INTRODUCTION 


Several researchers have examined the 
aeroelastlc stability of a helicopter blade 
In hover and forward flight (see recent reviews 
l" 3 ). The phenomenon is complex Involving non- 
linear structural. Inertial and aerodynamic 
forces. With a forward flight, the equations 
of blade motion get more Involved because of 
the presence of many periodic terms. Due to 
the complexities of formulation and analysts 
of rotorary-wlng dynamics problems, most of 
the analytical studies are of limited scope; 
more so, In forward flight conditions. The 
objective of the present paper Is to examine 
aeroelastlc stability In forward flight, 
Including the effects of dynamic Inflow on 
stability results. For this a simple flap- 
lag-torsion blade model consisting of three 
degrees of motion will be studied. 


For design and analysis of a helicopter 
rotor. It Is essential to analyze Its 
aeroelastlc stability. For this, a study on 
the dynamics of a single blade forms an 
important fundamental step to the complete 
understanding of the rotor-body dynamics. The 
blade stability analysis consists of three 
major phases; vehicle trim, blade steady 
response and stability of perturbation 
motion. The vehicle trim solution determines 
control settings for persrrlbed flight con- 
ditions and Is calculated from the vehicle 
overall equilibrium equations. The blade 
response solution consists of time dependent 
blade position and Is calculated fron the 
blade equilibrium equations. In the calcula- 
tions of blade response one needs the vehicle 
trim solution. For stability solution, a 
perturbation Is given to the blade at Its 
equilibrium position and the suKequerit 
response amplitude Is Investigated fur stabi- 
lity. For stability calculations, one needs 
tne vehicle trim solution as well as blade 
response solution. These three phases of 
study are Inherently coupled. A complete 
coupled solution Is very Involved and there- 
fore most of the researchers uncouple these 
three phases and study each phase separately. 
It Is possible however to achieve a certain 
degree of coupling between three phases 
through an Iterative process. 

The simplest form of a rotor blade 
representation Is the rigid blade model with 
spring restrained hinges. Many researchers 
have examined the aeroelastlc stability of 
this simple blade configuration. For 
example, Peter* and Kaza and Kvaternlk 5 
Investigated the aeroelastlc stability of 
two-degree flap-lag blade In forward flight. 
An Improvement for this type of modelling is 
to Introduce a third degree of motion, l.e. , 
feather rotation. A better representation 
for a hingeless blade Is to treat It as an 
elastic beam. As an example, Friedmann and 
Kottapalll 6 have Investigated aeroelastlc 
stability of flap bending, lag bending and 
torsion of an elastic hingeless blade In for- 
ward flight. In the present paper, a simple 
blade representation consisting of three 
degrees of motion, flap, lag and feather 
rotations, Is used to study the stability 
phenomena In forward flight. 

There are many forms of vehicle trim 
solutions available In literature. Johnson 7 , 
for example, presented In a summary form 
many trim options. For free flight con- 
ditions, the control settings and the 
vehicle angles are determined from the 
satisfaction of three force and three moment 
equilibrium equations. One of the popular 
trim procedure Is to neglect altogether 
yawing moment equilibrium equation and 
thereby neglect the Influence of tall rotor 
on solution. This form of trim solution Is 
used In the present paper. The next simple 
form of trim procedure® Is to neglect the 
lateral force equilibrium equation, and 
thereby exclude the determination of lateral 
shaft tilt angle (+.) from equilibrium 
equations. Generally, this may cause 
only slight Influence on trim and 


stability solutions, because the shaft 
lateral tilt angle does not introduce any 
vertical flow component on the blade. Some 
researchers 4 *®*® have simplified the trim 
procedure further by assuming that the 
vehicle center of gravity lies at the rotor 
hub, and thereby neglect the equilibrium of 
pitching and rolling moments of the vehicle. 

This will cause cyclic flap angles Bi c and 
Bis (# 1th aspect to h o plane) to be zero. 

Here, the control s- logs and shaft angle 
ot s are calculated from the vertical and 
longitudinal force equilibrium equations. 

This may again have a smail Influence on trim 
solution for free level flight conditions at 
low forward speeds. At high forward speeds, 
the cyclic flap angles are not small and 
therefore must not be neglected. Another 
form of trim procedure called moment trim Is 
often used by many researchers 4 *9* 10 , and for 
this the solution is calculated from the 
rolling moment and the pitching moment 
equilibrium equations. The force equilibrium 
equations are not considered. Here, the 
rotor cyclic controls (ai c and 9i s ) are 
calculated for a prescribed shaft angle 0 $, 

Some people refer it as a wind tunnel trim 
and It can be quite different from propulsive 
trim 4 * 6 . 

The blade time dependent poslton Is 
calculated from blade equilibrium equations. 
These are coupled equations and contain 
nonlinear geometric terms as well as periodic 
terms. The objective Is to calculate steady 
periodic response solution. In the present 
paper, the nonlinear equations are solved In 
the rotating frame In an Iterative procedure 
based on Floquet theory 11 . a somewhat similar 
type of quaslllnearlzatlon procedure was used 
by Friedmann and Kottapalll®, The solution 
contains all harmonics for flap, lag and tor- 
sion response amplitudes. Another popular 
method, harmonic balancing^ (Fourier Series) 

Is quite commonly used to calculate the blade 
steady response where response Is assumed 
periodic and consists of sum of finite har- 
monics. This procedure gets quite involved 
for coupled systems with nonllnearltles. 

Qul.e frequently, researchers 4 have obtained 
simple response solution using harmonic 
balance method where the flap response is 
assumed to undergo a single harmonic motion 
(So, Bi c and Sj s ) and the lag and torsion 
responses are ne lected. In literature 1 *®, 
the Importance of accurate determination of 
blade equilibrium position on blade stability 
has been pointed out. Including noi,1 1 r ar 
terms as well as higher blade harmonics. 

For stability analysis, t urbatlon 

equations of motion are llnaa’ oout the 
blade equilibrium position ana u.-se equations 
contain many periodic terms. These linearized 
equations are solvad using three different 
approaches In the present paper. The first 
approach® is to analyze the stability of the 
blade In the rotating reference frame using 
Floquet transition matrix theory. This approach 
Is applicable If the Inflow Is assumed to be 
steady. The second and third approach analyze 
the stability or rotor perturbation equations In 


the fixed reference frame. It Is assumed that 
the rotor Is tracked and all the blades are Iden- 
tical, The blade equations In the rotating 
reference frame are transformed to the fixed 
reference frame as rotor eguatlons using Fourier 
coordinate transformation^. In the present 
paper, these transformations are performed numeri- 
cally and thus the working through the laborious 
algebraic expressions Is avoided. In the second 
approach, the rotor equations In the fixed frame 
are solved using Floquet transition matrix theory. 
Through the coord, ate transformation, many 
periodic terms present In the rotating frame get 
cancelled out In the fixed frame. Therefore, the 
rotor equations in the fixed frame contain only 
selected periodic terms, for example, third 
harmonic for tnree-bladed rotor and second and 
fourth harmonic for four-bladed rotor. In the 
third aoproach 14 , a constant coefficient approxi- 
mation Is made by averaglrg out periodic ten s 
and solving the resulting equations. 

In all these three approaches an eigen- 
analysis is made and the nature of eigenvalues 
explains the stability of the blade. Another com- 
monly used method*® i s numerical Integration of 
complete equations. This approach is though 
simple in Implementation, but Is quite heavy from 
computation point of view. 

For trim and response solutions, the quasi- 
steady approximation Is used for the determination 
of aerodynamic loads. For the perturbation 
solution, the unsteady aerodynamics effects ran 
be Important and these are Introduced In an 
approximate manner, through a dynamic Inflow 
modelling. The effect of dynamic Inflow on 
coupled flap-lag two degrees :f motion In forward 
flight has been Investigated earlier’ -10 and has 
been shown f be quite Impo. tant for bl ac* stabi- 
lity. In the present paper, the influence of 
Inflow dynamics has been Investigated for a 
coupled flap-lag-torsion motion with improved trim 
and response solutions. The dynamic Inflow 
modelling Is based on the actuator disk tneory. 
This necessitates the transformation of blade 
aerodynamic forces to the fixed reference frame 
and therefore only second and third approach can 
be conveniently used to analyze blade stability. 

In the paper, the effects of several para- 
meters on bl de stability Is examined, including, 
structural cijpllng, pitch-flap coupling, pitch- 
lag coupling, lag stiffness, torsion stiffness, 
steady Inflow distribution, dynamic Inflow, blade 
response and constant coefficient approximation. 

EQUATIONS OF MOTION 

The blade Is assumed to undergo three degrees 
of motion: rigid body flap, lag and feather rota- 

tions about hinges at the t'adp lot, with hinge 
springs to obtain arbitrary natural frequencies. 
The hinge sequence is fl»p Inboara, lag, ano then 
feather outboard. The flap angle s Is positive up 
the lag angle c Is positive aft (opposite to 
rotation) and the feather angle e Is positive nose 
up. The equations of motion ere derived for this 
configuration, and in general, terms up to second 
order are retained In the flap and lap equations 
and terms up to third order a>’e retained In the 
feather equation. The equations are 
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Feather Equation: 
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where 1^ Is the ratio of the feathering inertia to 
flapping inertia; X, is the chordwise offset of 
the center of gravity from the pitch axis 
(positive aft), the cl and ?e are the viscous 
damping coefficients; the u g , u and u g are the 
nonrotating natural frequencies 5 of the blade 
(divided by rotational speed n) and 8 is the pre- 
cone angle. The Ma, M c , M e are f >.e aerodynamic 
flap, lag and feather moments r-s^ectively. For a 
uniform blade the nondimensionai rotating frequen- 
cies are given as 
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Kg, K. are the combined hub and blade stiffnesses 
in f’ap and lag directions. The R s *0 represents 
the configuration with blade part as rigid and all 
the flexibility concentrated at the hub. The 
structural fully coupled is represented by R s *l 
and this idealizes flexible blade with rigid hub. 
The intermediate values of R s represent the case 
where both blade as well as hub are flexible. 

Quasisteady airfoil characteristics are used 
tw obtain the aerodynamic forces. The pertur- 
bation section aerodynamic forces and pitch moment 
(in the shaft axis) are 
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where is the chordwise offset of the aerodynamic 
center from elastic axis (positive aft), c Is the 
chord, V is the resultant velocity and U p and 
are airflow velocity '■omponents in tangential and 
normal directions (Fig. 1(b)). The steady and per- 
turbation flow components for forward flight are 

Steady: 

* 

U t a fx(l-e) - x(l-e); + p sin ^ - cu cos 

Up - (X - x(l-e) 8c + x(l-e)B + uB cos 

Perturbation: (3) 

8U t » { 6 c’> cos ♦ - x(l-e)«c} 


The e Is the hinge offset (divided by radius of 
the blade) and R $ is «-h 2 structural coupling para- 
meter. A simple means of representing structure 
coupling effect in the rigid blade representation 
Is Illustrated In Fig. 1(a). It is used to charac- 
terize the hub to biade stiffness and Is defined as 
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iUp « [-x(l-e)s«c - x(l-e)c«S + x(l-e)«B + 

u«8 cos tp - - (! + ~) «e} 
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The perturbation aerodynamic moments required for 
the stability analysis a<e written as 

B 

* / x * 5F ^ • dx 


The perturbation inertia matrix M r P, damping 
matrix C r P, stiffness matrix JC r P also contain 
periodic terms. To determined lade stability one 
needs blade response solution. To calculate this 
blade response solution one needs the vehicle trim 
solution. 
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where is the nonci rculatory aerodynamic pitch 


VEHICLE TRIM SOLUTION 


The propulsion trim which is described here 
simulates the free flight condition. The trim 
solution in forward flight involves the calcula- 
tions of pilot-control setting as well as the 
vehicle orientation for a prescribed flight con- 
ditions. For a specified weight coefficient 
C and a fixed forward speed (u) the trim solution 


moment and is expressed as 
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(5) 
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where x is the nondimensional distance from the 
hub and u is the advance ratio (V cos a s /(fl). 

The final blade equations of motion in for- 
ward flight can be written as 

M r U)X* + C r (*)X r + K r ( *)X r = F NL (*,X r ,X r ) (6) 


where the inertia matrix M r< the damping matrix 
£ r and stiffness matrix K r contain periodic terms. 
The vector X r consists o'? three states; flap, lag 

and torsion deflections in the rotating system. 

The (*) shown in the equations refers differen- 
tiations with respect to * and ip is the azimuth 
angle (nondimensional time, (it). All the 
geometric nonlinearities are put into the force 
vector Jnl, The blade response is calculated 
from the solution of above equations. 

For the stability solution, the flutter 
motion is assumed to be a small perturbation 
about the blade equilibrium position. 


X. ■ X + 
— r ~r 


~r 


( 7 ) 


The final linearized perturbed equations are 
obtained as 


M r P (*)« r K r P( Mr'ir ) V^^'ir'Jr^r ‘ 0 (*> 


evaluates fy, 8 ls . %, \ c . ® 1$ , <^ p . * s and 

X. The trim solution is calculated from the 
vehicle equilibrium equations. F’g. 1(c) shows 
the forces and moments acting on the vehicle. 

Vertical force equilibrium: 


W-T cos (a - 8pp) cos 4 S + D sin 9pp - H sin 


(a - 6j. p ) + Y sin * s = 0 

(9) 

Longitudinal force equilibrium: 


D COS t^p+H cos (o-fl^p)-T sin (a - 6p p ) * 

0 (10) 

Lateral force equilibrium: 

Y f + Y cos * s + T sin = 0 

Pitching moment: 

HD 

M y + M yf ' w ^ X cg - h sin “) * D cos 

(a + ^ p )h - D sin (a + e pp) *X cg =U 

(12) 

Rolling moment: 


M x + M xf * V h cos ♦* + Y f Y cg Sin ♦* 


+ Wh sin ♦j - W Y cg cos ^ = 0 

(13) 

Where T is the rotor thrust and Y and H are the 
side and drag forces. These are five vehicle 
equilibrium equations. For trim solution one 
also needs rotor equilibrium equations in 
simplified form. These are 

f 2 * (Flapping equation) dip = 0 
2* 0 

(14) 

-i- j 2 " (Flapping equations) cos tp d ip * 0 
Z* 0 

(lb) 

-i- j 2 * (Flapping equations) sin ip d ip ■ 0 

2* 0 

116) 


2J(5 


Tne induced inflow Is related to the rotor thrust 
as 

Ct 

> * u tan a + (l+k *x cos <r*-K «x sin *) 

2/(i 2 +x 2 * y ( 17) 

where K x and are obtained from Drees 
model ^ an( j <(5 expressed as 

K x > 4/3[(1-1.8 u 2 )/1 - ( Vu) 2 - V»] 



For hover K * K =0 

x y 

These are nine equations with nine unknowns ( Bjy. 

B lc* hs* V 0 lc* °ls* X * °’V and these are 
solved numerically by Iterative method. 

BLADE RESPONSE SOLUTION 

The blade response solution involves the 
determination of time dependent blade position and 
is calculated from the blade equations (6). For 
the calculation of response solution one needs 
vehicle rim solution. The nonlinear response 
solution is obtained from nonlinear periodic 
equations (6) using Floquet theory 11 . These 
equations are expressed in the state vector form 

j r * A(*)I r * G r (*. Y r , Y r ) (18) 

where Y r is the state variable vector involving 
six states. 

First a linear solution is calculated after 
dropping all nonlinear terms. For this the 
initial conditions are calculated from Floquet 
theory as 

Y(0) « (I_ - J3(2w))- X r e (2«) (19) 

where Y e (2t) is the complete solution after one 
revolution with rest initial conditions and y(2v) 
is the Floquet transition matrix. For numerical 
Integration of the equations, a fourth order Runge- 
Kutta algorithm is used. The next step is to 
obtain the Initial conditions for the complete 
nonlinear problem. This is done in an iterative 
manner. As a first guess, the above linear solu- 
tion is used as an initial ..ctor Y.(0) for the 
nonlinear solution and the complete* response 
Y e (2») after one revolution is calculated. The 
updated Floquet transition matrix £ Is function of 
response amplitude and Is calculated by perturbing 
the estimated Initial conditions _Y (0) by a smail 
perturbation vector e. 

Q(Zx) « C— a (1) (2x) - Y e (2*)), I (yJ 2) (2v) 

e l h 

- T (2i.)),.,,-i(Yj 6) (2«)-Y (2v))3 (20) 

6 e 6 

Where Y^(2») is the response with Initial con- 
ditions of 



So a new set of Initial conditions for the nonli- 
near analysis are obtained 

T(0) * Y e (0) + [I - Q]' 1 (Y e (2«) - Y e (0)) (21) 

This procedure is repeated till a converged set of 
initial conditions is obtained. Typically it 
takes about 2 to 3 iterations to obtain converged 
solution. Once the initial conditions X(°) are 
obtained, then the total response Y(») for any 
time in a revolution is calculated numerically 
using time integration (Runge Kutta). This gives 
us the nonlinear equilibrium deflection of the 
blade along the azimuth. 

FLUTTER SOLUTION 


The linearized perturbed equations of motion 
(8) is written in the state vector form as 

s !r ‘ A( *r* Jr* *^r (22) 


where Y r , J r are the blade equil librium position 

in the rotating reference frame and «Y and 4Y 
are perturbation states. These linearized r 
equations are solved for stability using Floquet 
transition matrix theory 1 ^. Here the eigenvalues 
of transition matrix of A c-n be written in 
characteristic exponents -•m 

\ * \ + n (23) 

and the mode in stable when a. < 0, 

For the perturbation solution the unsteady 
aerodynamics effect can be important and these are 
introduced in an approximate manner, through a 
dynamic Inflow modelling. The wake inflow is 
perturbed about the steady inflow X 


X » X + «x (24) 

where 6X is the perturbed inflow 
component. 

A linear variation of perturbed inflow 
is used 

5x • Sky + 6Xj C x cos * + «Xj s x sin <> (25) 

The dynamic Inflow components «x 0 , 

4 *1C> 4x ls are related to rotor unsteady 
aerodynamic forces and moments 
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The «C T , aC^ , ac ffl are the perturDed 
x y 

thrust, roll moment and pitch moment and 
these are obtained for the 1th olade. 

The m and t-1 matrices used here are 
evaluated analytically based on the 
actuator disk theory in Ref. 16. The 
nonzero elements of m and t are 


m n » —— m~~ * m,, - - ■i£- 

H 75, cl 33 45, 

, , 1 , . *33 = _ 4 l 

*11 77 *22 sfna TFsTna v 
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u 2 +X(X+Ai) 

/u 2 +X 2 


and 
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tan 


-1 
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1 

— 5° (Differential equation) 


If m=l 
b 



These transformations are carried out 
numerically. The final equations of 
motion In fixed co-ordinate frame can be 
written as 


«Y p » A f (*) + N f (*)4X (29) 


The unsteady force «F In dynamic Inflow 
equations Is given as 


€ -li(«fi F + 5 2 (*)“ (30) 

Putting together the rotor equations 
with the inflow equation, one gets 



“ 



Nf (*) 


A 

M'V-M'V 1 



- 




(31) 


X, Is the induced inflow due to steady 
rotor thrust. This model gives a quite 
accurate description of dynamic inflow 
as concluded in reference 3. 

The disk loading Is approximated In 
terms of the blade loading, 6F Z as 

~ K 1 

* T * — l f dx 

yN b k-1 e 2 K 


The above equations in the fixed frame 
contain only selected harmonics, for 
example third harmonic for 3 bladed 
rotor, second and fourth for 4 bladed 
rotor. These linearized periodic 
equations are solved using Floquet tran- 
sition matrix theory and constant coef- 
ficient approximation approach. In the 
constant coefficient approximation 
approach the periodic terms are averaged 
out by applying the operator 



N b 1 

— i / (^a * dx s1n \ 

yN.k-1 e X X 

b (27) 
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■7- / ( ) d* (32) 
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«C M » - — l I (#■ ' ). xdx cos * 
"y yN b k*l e 2 k ^ 


With the Inclusion of dynamic Inflow, It 
is convenient to analyze blade stability 
In the fixed reference frame. The 
coupled blade equations are transformed 
to the fixed reference frame from 
rotating frame using Fourier co-ordinate 
transformation! 3 . For f 0ur bladed rotor 

<v 4 > 


and then solved as an eigenvalue problem. 


RESULTS AND DISCUSSION 


The aeroelastic stability is examined for 
a four-bladed rotor with Lock number y * 5.0, 
solidity ratio o ■ 0.05, feather inertia to 
* 

flap ratio I { »0,0003 and with zero precone. 

The blade offsets such as the Xj , the chord- 
wise center of g vity offset from pitch axis 
In terms of radius and X 4 aerodynamic centre 
offset from elastic axis are set to zero. The 
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fuselage centre of gravity lies on the shaft 
axis and Is assumed to be at a distance 0.2R 
below the rotor center. The aerofoil charac- 
teristics used are 


V 5 - 7a 

C d = 0.01 

C = -0.02 

m 

The helicopter drag coefficient In terms of 
flat plate area ratio (f/*R 2 ) of .01 is used. 
The blade flap and torsion frequencies 
(rotating) are 1.15/rev and 5.0/rev respec- 
tively. Two different lag frequencies are 
used; 0.57/rev for soft Inplane rotor and 
1.4/rev for stiff inplane rotor. The soft 
inplane rotor conf igurat ion was taken as a 
matched stiffness case (« s * u^). 

First, results were calculated for some 
selected cases to make comparison with those 
cf ether authors for Identical conditions. 

Tha vehicle trim was calculated for uniform 
inflow conditon with the center of gravity 
lying at the rotor hub (h=0) and the results 
obtained were identical to those of Ref. 18. 
The blade stability of two-degree flap- lag 
blade was calculated using simple response 
solution (single flap harmonic) and for uni- 
form steady inflow condition. The lag damping 
values calculated for various flight con- 
ditions were quite Identical to those of Ref. 
4. The inclusion of dynamic inflow on the 
stability of this two-degree motion blade was 
checked with the results of Ref. 9, and again, 
the comparison was quite satisfactory. Then, 
the flap-lag-torsion stability results were 
calculated for soft and stiff Inplane con- 
figurations for steady Inflow conditions. 

For perturbation solution, a nonlinear 
equilibrium position was used. The stability 
results showed some comparable trends with 
those of Ref. 6, in which results are obtained 
with impioved structural modelling for the 
blade (elastic beam). 

Figs. 2(a) and 2(b) show the vehicle trim 

solutions for of .2 and .1 respectively. 

o 

The propulsive trim parameters So, eje, 8i s , 
“HP. *s. ®0> BlC • Sis are plotted for dif- 

ferent forward speeds (In terms of advance 
ratio u=V/fjp) . The solution is calculated 
iteratively from nonlinear equilibrium 
equations (large angles). The flight path 
angle e™ is assumed to be zero. For steady 
inflow, a linear distribution model (Drees) Is 
used. These trim parameters are defined in 
the hub plane axes system and so the cyclic 
flap angle 8, c is small for even large u. For 
larger y, theShaft has to tilt more to com- 
pensate the Increase In parasite drag and 
hence an,, increases with y and thereby cause 
Inflow x to Increase. For large thrust levels 
Cy/n, control requirements are large but the 
angle o HP is small. The reason for lower 
shaft angle at high C w / 0 is due to the fact 


that, for same parasite drag (for a given for 
ward speed) the shaft has to tilt less to 
balance the parasite drag as the thrust magni- 
tude is more. It is also observed that trim 
calculated from linearized equilibrium 
equations (small angles assumption) is quite 
close to the above nonlinear solution except 
for large advance ratios. The Influence of 
Inflow distribution is primarily on longitudi- 
nal and lateral cyclic pitch e^j, e^. The 

Cw 

effect of uniform inflow for — = 0.1 is shown 

o 

by dotted line in Fig. 2(b). The effect of 
inflow distribution on other trim parameters 
is quite small. 

Figs. 3(a)-(c) present time dependent 
equilibrium position of blade for one 
complete revolution. The response solution 
in terms of flap, lag and torsion deflections 
(angles) is calculated iteratively from nonli- 
near blade equations (rotating frame) using 
the floquet theory described earlier. These 

. C w 

results correspond to a — =■ 0.1 and advance 

ratio y = 0.2. For comparison, the linear 
response solution (dotted) is also presented. 
The geometric nonlinearlties are important for 
lag response and play less important role in 
flap and torsion response calculations. 

For numerical results, a convergence 
study was conducted to determine time steps 
needed in one revolution for time integration 
(Runge-Kutta) for both response as well as 
stability calculations using Floquet Theory. 

It was concluded that 120 time steps are quite 
adequate for well converged (four 
significant-digit) response and stability 
solutions. For stability results, only the 
lowest damped lag mode is presented. The flap 
and torsion modes are comparitively high 
damped modes and are not presented here. The 
damping In terms of real part of complex 
eigenvalue, aq Is shown. Note a, » CL“c> 
where damping ratio of lag mode and w,/g Is 
frequency of lag mode nondlmenslonall zed with 
respect to rotational speed. 

In Fig. 4, the effect of torsion 
flexibility on the blade stability is 
shown. For these solutions, the simple 
blade response Is used and the stability 
roots are obtained In the rotating frame 
with steady uniform inflow conditions. 

There is a disparity between two results, 
clearly showing the Importance of Inclusion 
of torsion flexibility for blade stability 
analysis. This has been pointed out by 
other authors 6 . 

Fig. 5 shows the effect of blade 
response solutions on statillty. The lag 
mode damping is calculated using three 
types of blade equilibrium solutions; these 
are complete nonlinear solution, linear 
solution and simple solution. The linear 
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and nonlinear solutions contain all har- 
monics for flap, lag and torsion modes 
whereas the simple solution consists of 
single flap harmonic only. These solutions 
are obtained In the rotating frame with 
steady uniform Inflow conditions. Though it 
Is computationally less Involved to use 
simple response solution, the results are 
poor In accuracy. For accurate results It 
Is needed to use a complete nonlinear blade 
response solution. 

Fig. 6 presents the Influence of 
steady Inflow distribution on the blade 
stability. These solutions are obtained In 
the fixed frame with dynamic Inflow Included. 
In the figure the damping of the lowest 
damped, low frequency cyclic lag mode 
(regressive mode) Is presented. Two types 
if steady Inflow model are used; uniform 
distribution and linear distribution 
(Drees). The uniform distribution 
underpredicts lag damping. For subsequent 
results, the linear Inflow model (Drees) is 
used. 

Figs. 7 and 8 present the damping of 
low frequency cyclic lag mode for different 
advance ratio u. Three sets of results are 
shown and these respectively represent 
dynamic Inflow Floquet results (full line), 
dynamic Inflow constant coefficient 
approximation results (big dots) and steady 
Inflow results (small dots). In Figures 
7(a) and 7(b), the stability results are 
shown for stiff Inplane rotor for C w / 0 - .1 
and .2 respectively. For this case the low 
frequency mode Is a regressive mode. The 
constant coefficient approximation Is quite 
satisfactory for small advance ratios u. 

The Inclusion of dynamic Inflow Is Impor- 
tant for low forward speeds. This shows 
that for large forward speeds unsteady 
aerodynamic effects are not important and 
quasisteady approximation Is quite adequate 
for blade stability analysis. At large 
thrust levels the Influence of dynamic 
Inflow and constant coefficient approxima- 
tion Is large because of larger aerodynamic 
forces involved. Figs. 8(a) and 8(b) 
show the blade stability results for soft 
Inplane rotor for C w / 0 of .1 and .2 respec- 
tively. For this case the low frequency 
mode Is a progressive mode. Again for the 
matched stiffness configuration, the 
results are quite Identical to the stiff 
configuration. Constant coefficient 
approximation less satisfactory for high 
advance ratios, the dynamic Inflow incl - 
slon Is more Important for low advance 
ratios. 

In Figs. 9 and 10 the effect of struc- 
tural coupling on lag mode stability is 
presented for stiff Inplane and soft 
Inplane rotors respectively for Cw/ 0 of 
0.1. The earlier results were calculated 
for blades with no structural coupling 
(R««0). This Idealizes the configurations 
with blade part as rigid and all the flexi- 
bility concentrated at the hub. The struc- 
turally fully coupled Is represented by 


R«*l and this idealizes flexible blade with 
rigid hub. The value of R s less than 1 
represents Intermediate cases where both 
blade as well as hub are flexible. The 
stability results are obtained using 
Floquet theory with dynamic inflow effects 
Included. The structural coupling has an 
Important effect on blade stability for 
stiff inplane rotor. In fact with large 
structural coupling the blade becomes more 
stable. This is because with a large 
structural coupling R , the weakly damped 
lag mode gets coupled s with well damped flap 
mode and thereby stabilizes the lag mode. 
The effect of structural coupling is negli- 
gible on the soft lag configuration and 
this is because the configuration con- 
sidered is matched stiffness case. On 
these figures the results are also plotted 
with steady linear inflow aerodynamics. 

For stiff lag rotors, the effect of inclu- 
sion of dynamic inflow is large for con- 
figurations with zero structural coupling, 
and the Influence Is quite stabilizing. 

With large structural coupling R the 
effect of dynamic inflow is less s and it is 
destabilizing. 

Figs. 11 and 12 show the effect of 
torsional stiffness on lag mode stability 
for stiff Inplane and soft inplane rotors 
for C w / p of .1. Results are obtained for 
three different torsional frequencies and 
these are 2.5, 5 and 10 per revolution. 

For both soft inplane and stiff inplane 
rotors, Increasing torlonal stiffness 
Increases lag damping (more stable) for 
lower forward speeds (low u) and decrease 
lag damping (less stable) at higher forward 
speeds. 

The effect of pitch-flap and pitch-lag 
coupling terms on blade stability is 
studied by modifying the feather angle in 
the flap-lag equations 
(two-degree-of-freedom) 


Vf ■ 6 - \ « - \ ‘ 

The pitch-flap V.vvip • irq Kp. is positive 
flap up/pitch down, and tho pitch-lag 
coupling is positive lag back/pitch down. 
These couplings are caused due to torsion 
dynamics or kinematic couplings. Figs. 
13(a) rod 13(b) show the influence of 
pitch- ug coupling on lag mode stability 
for stiff Inplane and soft Inplane rotors 
respectively for C w /o of 0.1. The positive 
pitch-lag coupling stabilizes the low fre- 
guency cyclic lag ■"ode for stiff inplane 
rotors, and destablizes this lag mode for 
soft Inplane rotors. The opposite effect 
Is seen with the negative pitch-lag 
coupling. A similar type of observation Is 
made for hovering blade stability in Ref. 
17, Figs. 14(a) and 14(b) show the effect 
of pitch-flap coupling on lag mode stabi- 
lity for stiff Inplane and soft Inplane 
rotors respectively for C w /o of 0.1. A 
negative pitch-flap coupling reduces the 
flap frequency, and it produces a 
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stabilizing effect on lag mode for low for- 
ward speeds and a destabilizing effect at 
higher forward speeds. A positive pitch- 
flap coupling raises the flap frequency, 
and it has a comparatively small effect on 
lag mode stability. Also, It can be seen 
the effect of pitch-flap coupling on blade 
stability is much smaller as compared to 
that of pitch-lag coupling. 


CONCLUSIONS 

An aeroelastic stability of a simple 
three-degree-of-freedom blade model In for- 
ward flight is examined. The nonlinear 
time dependent blade equilibrium position 
is calculated using a quasi linearization 
procedure based on Floquet theory. The 
perturbation solution is obtained using 
Floquet transition matrix theory as well as 
constant coefficient approximation in the 
fixed reference frame. The stability 
results are calculated for both stiff- 
inplane and soft-inplane blade con- 
figurations. The inclusion of torsion 
degree of motion is important for blade 
stability. The nonlinear time dependent 
periodic blade response has a significant 
influence on blade stability. For steady 
inflow distribution, the linear variation 
(Drees) is somewhat stabilizing for lag 
mode damping as compared to uniform distri- 
bution. The effect of dynamic inflow on 
lag mode stability Is small at high forward 
speeds (u>.3). The constant coefficient 
approximation appears satisfactory for low 
forward speeds ( y < .2) . The structural 
coupling produces stabilizing effect on 
blade stability for stiff lag rotors. For 
matched-stiffness configurations, there is 
no effect of structural coupling on blade 
stability. Raising of torsional stiffness 
Increases lag mode damping (more stable) at 
lower forward speeds ( u < .15) and decreases 
lag damping (less stable) at higher forward 
speeds (u > .15). A positive pitch-lag 
coupling stabilizes the low frequency 
cyclic lag mode for stiff lag rotors and 
destabilizes this mode for soft-lag rotors. 
An opposite effect is seen with negative 
pitch-lag coupling. The effect of pitch- 
flap coupling on lag mode damning Is small 
as compared to pitch-lag coupling effect. 
The negative pitch-flap coupling stabilizes 
lag mode at low forward speeds and destabi- 
lizes it at higher forward speeds. The 
positive pitch-flap coupling has a little 
influence on lag mode stability. 
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Fig. 1(a) Spring model for elastic blade and hub 



Fig. 1(b) Blade section Aerodynamics 
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Fig. 1(c) Helicopter In forward flight 

showing vehicle trim configuration 
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Fig. 3(c) Torsion equilibrium angle 


Fig. 4 Effect of torsion degree of freedom on 
lag mode stability for C w /o=0.2 
(v =1.15, v =1.4, R =0, steady uniform 

P g 5 

inflow) 
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Fig. 5 Effect of blade response solutions on 
lag mode stability for C w /o*0.2 
(v *1.15, v *1.4, v *5, R *0, steady 

P g 0 5 

uniform inflow) 


Fig. 6 Effect of steady Inflow on low frequency 
cyclic lag mode (Regressive) stability 
for C w /o*0.1 (v e -l .15, v c *1.4, v e *5, 

R *0, Dynamic Inflow) 



Fig. 7(a) Stiff-lag rotor, C w /o=0.1 Fig. 7(b) Stiff-lag rotor, C w /o=0.2 

Fig. 7 Damping of low frequency cyclic lag mode for different advance ratio 
y for a stiff lag rotor (v =1.15, v =1.4, v =5, R =0) 

C 0 5 



Fig. 8(a) Soft lag rotor, C^/a^O.l Fig. 8(b) Soft lag rotor, C w /o=0.2 

Fig. 8 Damping of low frequency cyclic lag mode for different advance ratios 
of soft-lag rotor (v 8 *1.5, v ? »0.57, v g *5, R s »0) 
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Fig. 9 Effect of structural coupling on low 
frequency cyclic mode (regressive) for 
stiff Inplane rotor for C w /o*.l 
(v £ =l .5, v =1 .4, y 0 =5) 



Fig. 11 Effect of torsional stiffness on low 
frequency cyclic mode for stiff lag 
rotor (C^/o*.!, v e «1.15, v «1.4, R s *0, 
Dynamic Inflow) 


Fig. 10 Effect of structural coupling on low 

frequency cyclic mode for soft lag rotor 
(C w /a=.l, v 0 =*l .15, v c =0.57, v 0 =5) 
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Fig. 12 Effect of torsional «Hf f ness on low 

frequency cyclic mode for sort lag rotor 
(C t ,/o“.l, v *1.15, v *.57, R »0, Dynamic 

r. p 

Inflow) 
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DISCUSSION 
Paper No. 16 

FLAP-LAG-TORSION STAB I L ITT IN FORWARD FLIGHT 
Brahmananda Panda 
and 

Inderjlt Chopra 


Wayne ■ ’ nnson, NASA Ames Research Center : Could you describe again the manner in which you 

trinmed che rotor in forward flight? I noticed in your plot of flapping motion that you had 
about one degree of 1 per rev flapping motion, so could you tell us how the rotor was trimmed 
for these results? 

Panda : First, we obtained the vehicle trim equations and calculate the control 9 0 , 0 lc , and 

9 1s . And using that with the coupled flap- lag- tors ion blade equations for solving t.ie trim. 

Johnson: Did you trim to zero moment about some center of gravity of a helicopter below the 

rotor hub, is chat what you did? 

Panda: Yes. 

Peretz Friedmann, University of California. Los Angeles : I wanted to be sure that I understand 

what kind of model you are using. You have an offset-hinged, spring-restrained blade model with 
flap, lag, and torsion degrees of freedom, is that correct? 

Panda : Yes. 

Friedmann : Then I Just wanted to comment that at the European Rotorcraft Forum this August a 
gentleman by the name of Neelakanthan did exactly the same problem. He also did the elastic 
blade with two flap, two lag, and two torsional degrees of freedom and his results indicated 
that the model you are using is not safe at all times, so I am just suggesting that maybe you 
should qualify your conclusions. 

Bill Bousmai , U.S. Army Aeromechanics Laboratory : I haven't had time to look at these last two 

papers, but I think that most of the stuff Peretz has done over the years and other people — Dave 
Peters— ha-: not shown an 'notability at high speed and there is none shown here in the lag 
mode. It's Just the continuing stabilizing effect as the inflow Increases. Dr. Reddy showed 
Ploquet splitting roots and a destabilizing effect and I guess it's an open question [whether] 
this is something new. Was it Just that his torsion [frequency] was 3 per r?v or what? 

Friedmann : I Just wanted to say. Bill, that the instability that Reddy has suown is one which 

appeared for a stiff inplane case in the paper which kottapalli and I have written and he Just 
checked it out and he got t.ie same instability. 


258 


